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A new calculation of the mass fraction of primordial black holes 
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We revisit the calculation of the abundance of primordial black holes (PBHs) formed from primor- 
dial density perturbations, using a formation criterion derived by Shibata and Sasaki which refers 
to a metric perturbation variable rather than the usual density contrast. We implement a deriva- 
tion of the PBH abundance which uses peaks theory, and compare it to the standard calculation 
based on a Press-Schechter-like approach. We find that the two are in reasonable agreement if the 
Press-Schechter threshold is in the range Ath — 0.3 to 0.5, but advocate use of the peaks theory 
expression which is based on a sounder theoretical footing. 



PACS numbers: 04.70.-s, 98.80.-k 



I. INTRODUCTION 



Primordial black holes (PBHs) may have formed dur- 
ing the early Universe, and if so can have observational 
implications at the present epoch, either from effects of 
their Hawking evaporation or from a contribution to the 
present dark matter density P, I3- That there is no un- 
ambiguous observational evidence of PBHs is a signifi- 
cant constraint on some possible types of early Universe 
physics. In particular, they are the only known way of 
constraining the density perturbation spectrum on ex- 
tremely short scales, and indeed until fairly recently pro- 
vided the most powerful upper limit on the spectral index 
of perturbations with an exactly power-law power spec- 
trum. 

However, the abundance of PBHs formed from a given 
initial power spectrum remains uncertain. The tradi- 
tional calculation takes the same form as the Press- 
Schechter calculation much used in large-scale structure 
studies 01 , where the density field is smoothed on a mass 
scale M (in this application taken to be at the time of 
horizon crossing), and those regions where the density 
contrast exceeds a threshold value Ath are assumed to 
form PBHs with mass greater than M. However the 
correct value for the threshold is quite uncertain. The 
'standard' value of 1/3 for a radiation-dominated Uni- 
verse was derived by Carr Q (see also Ref. Q), but was 
probably only ever intended as an order-of-magnitude es- 
timate. Subsequently, Niemeyer and Jedamzik ^ car- 
ried out numerical simulations of the collapse of isolated 
regions and found the threshold for PBH formation, in 
terms of the relative excess mass within the horizon, to 
be AM/Mh = 0.7. However, Shibata and Sasaki 
have pointed out that they formulate their initial data 
after horizon crossing, and hence their criterion cannot 
be related to the initial perturbations produced by, for 
instance, a period of inflation. 

More recently, Shibata and Sasaki devised a new 
approach to the formation of individual PBHs, seeking 
to find criteria on the metric perturbation rather than 
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the density field, and in a form which can be applied 
to superhorizon initial perturbations. They were able to 
specify a criterion in terms of whether the initial cen- 
tral value of a particular metric perturbation variable tp 
exceeds a threshold value. In this paper, we investigate 
the implications of this result for the abundance of PBHs 
formed. 



II. THE PBH FORMATION CRITERION 

We briefiy describe the PBH formation criterion of Shi- 
bata and Sasaki, whose paper can be consulted for the 
full details 6]. They define a metric variable ■0 from the 
spatial part of the metric on uniform-expansion hyper- 
surfaces [their Eq. (2.2)] as 



(1) 

where jij is the metric of the spatial 3-sections (through- 
out we assume a flat background and only consider scalar 
perturbations). Shibata and Sasaki numerically explored 
a range of initial configurations, all spherically symmet- 
ric, for the metric variable -0 in a radiation-dominated 
Universe, and were able to show that the central value of 
-0, denoted ipo, was a good indicator of PBH formation. 
They found that PBH formation took place provided V'o 
exceeded a threshold value V'o,th- The precise value of 
this threshold depended on the environment of the ini- 
tial configuration, and lay in the range from 1.4 for a 
density peak surrounded by a low-density region, to 1.8 
for a peak surrounded by a flat Friedmann-Robertson- 
Walker (FRW) region. 

We wish to relate the Shibata-Sasaki threshold crite- 
rion to quantities given in standard linear perturbation 
theory, where the spatial part of the metric tensor is given 
by J,] 



= [(1 + 27^) 5y + 2^^^jH^\ 



(2) 



where TZ is the curvature perturbation and iJx represents 
the anisotropic part. The gauge-invariant curvature per- 
turbation on uniform-density hypersurfaces C is defined 
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as 



p 



(3) 



with iJ, p, and <5p denoting the Hubble parameter, 
background density, and perturbed density, respectively, 
which then gives 



5y 
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20 



(4) 



on superhorizon scales, where the anisotropic part Ht 
is negligible. Note that on large scales uniform-density 
hypersurfaces coincide with uniform-expansion hypersur- 
faces (also known as uniform-Hubble hypersurfaces) 

A reasonable prescription for relating C and ■0 in the 
quasi-linear regime is 



exp(2C)=^^ 



(5) 



since by definition -0"* = exp (2A7V) where AA'' is the dif- 
ference in e-foldings between uniform-expansion hyper- 
surfaces, and we can argue that the uniform-expansion 
and uniform-density slices are almost equivalent even 
in the non-linear regime, so that C — AN 8, 9]. Us- 
ing Eq. 0, we find that the threshold values of ipo 
("00, th = 1-4 and 1.8) correspond to thresholds on ( of 
Cth = 0.7 and 1.2 respectively. 



III. THE PBH ABUNDANCE 

The observational constraints on the fraction of the 
energy density of the Universe in PBHs at the time they 
form, Qp-Bn{M) , can be very roughly summarized as 
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on any interesting mass scale. Detailed examination of 
particular constraints can give more accurate values for 
the limits at particular masses P, , but for our present 
purpose we need only have an approximate guideline. In 
any event, PBH formation calculations remain uncertain 
enough that high-accuracy observational constraints are 
unnecessary; nevertheless the production rate is normally 
so sensitive to quantities we might wish to constrain, such 
as the density perturbation amplitude, that useful con- 
straints can be extracted even from quite approximate 
calculations and constraints. 



A. Review of the standard calculation 

The traditional PBH abundance calculation 
(e.g. Refs. P, 01) refers to a quantity which in 
modern terminology would be known as the density 
contrast on the comoving (velocity-orthogonal) slic- 
ing, which we denote by A. The density contrast is 
smoothed on a scale R, and the calculation simply 



integrates the probability distribution P(A(i?)) over 
the range of perturbation sizes which form PBHs: 
Ath < A(i?) < Acut, where the upper limit arises since 
very large perturbations would correspond to a separate 
closed universe in the initial conditions ll j. In practice 
P(A(i?)) is such a rapidly decreasing function of A(i?) 
above Ath that the upper cut-off is not important. The 
threshold density is taken as Ath > w, where w = p/p is 
the equation of state 1]. This cannot of course be valid 
in the limit w ^ 0, but is thought to be acceptable for 
the radiation-dominated case w = 1/3 which is the main 
one of interest. 

The smoothed density contrast A(i?, x) is found by 
convolving the density contrast A(x) according to 



A{R, x) = V-^ J W{\^' - ^\/R)A{x')d^x' . 



(7) 



where R is the smoothing scale, W{y) is the window 
function used for the smoothing and V is the volume 
of the window function. If the initial perturbations are 
gaussian, this property will be inherited by the smoothed 
density perturbation so that 



P(A(i?)) = 
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where a^iR) is the variance of A(i?, x). 
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Here V^^k) = (A:'^/27r^)(|Afcp) is the power spectrum of 
A and W{kR) is the volume-normalized Fourier trans- 
form of the window function used to smooth A. It is not 
obvious what the correct smoothing function to use is; 
a top-hat smoothing function has often been used in the 
past 01 although it is sensitive to scales well within the 
horizon, which requires careful treatment 0| . We prefer 
to use a gaussian window function: 



W{kR) = exp 



(10) 



On comoving hypersurfaces there is a simple relation 
between the density perturbation and the curvature per- 
turbation (e.g. Ref. J. 3]) 



A{t,k) 
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where TZc is the curvature perturbation on comoving hy- 
persurfaces, which coincides with the curvature perturba- 
tion on uniform-density hypersurfaces, Eq. on large 
scales. The power spectra are related by 



VA{k,t) 



4(1 



(5-l-3w)2 



(12) 



3 



Then at horizon crossing we have 



(13) 



The fraction of the Universe which exceeds the thresh- 
old for PBH formation A{M) > Ath when smoothed on 
scale M, and hence will form a PBH with mass > M,^ is 
given as in Press-Schechter theory by 

/•OO 

^^PBH,Ps(Ath > M) = 2 / P(A(M))dA(M) 
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In this expression we have followed the usual Press- 
Schechter practice of multiplying by a factor 2, which 
can be thought of as allowing for the fact that the PBH 
formation happens in regions which are overdense with 
respect to the mean cosmological density. 

For the purpose of specific calculations in this pa- 
per, we assume a power-law primordial power spectrum 
VnAk) = AnAk/koY-^, so that 



al{M) = 



2{l + wY AnT[{n-l)/2] 



(15) 



(5 + 3w)2 [koRY-^ 

Spergel et al. [l| found, from the WMAPext-|-2dFGRS 
dataset, that An, = (0.8 ± 0.1) x 2.95 x lO^^ for = 
0.05Mpc"\ 

B. A new calculation using peaks theory 



The Shibata and Sasaki PBH formation criterion is ex- 
pressed in terms of the peak value of the fluctuation, ip^^ 
at i = (equivalently, at some early time when the per- 
turbation is on superhorizon scales, since tp is constant 
on superhorizon scales). Rather than Press-Schechter, it 
is therefore best suited to a calculation of the mass func- 
tion using the th eory of peaks, as extensively described 
by Bardeen et al. jla| ■ We will apply peaks theory to the 
initial value of the variable C- 

After smoothing the density field on a scale M, the 
number density of peaks with height greater than 
where v = C,t\i/ <yc,{M), is given (for high peaks) by [TslIT^ 



Speaks (i^, M) = 



(2^)2 V 3 ) 



l)e-V^ 
(16) 



^ Throughout we assume for simphcity that the PBH mass is equal 
to the horizon mass Afjj corresponding to the smoothing scale. 
This is not strictly true (and in fact the PBH mass appears to de- 
pend on the size and shape of the perturbations however, 
this uncertainty is not important when applying PBH abundance 
constraints, due to their relatively weak mass dependence. 



where {k'^){M) is the second moment, with respect to k, 
of the power spectrum 



{k'){M) 



(17) 



For a power-law power spectrum V({k) — A(^{k/ko)"~^ 
(with A^ = All, since on super-horizon scales C = T^c) 
and a gaussian window function, we have 



{k')(M) = 



n- 1 
2i?2 ■ 



(18) 



The number density of peaks with height greater than 
v, when smoothed with a gaussian filter on scale M , is 
then given by: 



where 



2(fcofl)"-i 
A^V[{n ~l)/2] 



1/2 



th ■ 
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The number density of peaks is related to the frac- 
tion of the Universe in peaks above the threshold by 
f^PBH, peaks (i^, > M) = n^caks{v,M)M/ p. Here M is the 
mass associated with the filter (which for a gaussian win- 
dow function is given by M ~ p{2'k)'^/'^R'^) so that 
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where 
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C. Comparison 



To use our results, we need to relate the comoving 
smoothing scale R to the horizon mass. The main case 
of interest is radiation domination, where w = 1/3. The 
horizon mass is given by 



An 



-1\3 



(23) 



when the scale enters the horizon, R — (aH) ^. During 
radiation domination aH oc a^^, and expansion at con- 
stant entropy gives p cx ^^^a~'^ ^| (where g^, is the 
number of relativistic degrees of freedom, and we have 
approximated the temperature and entropy degrees of 
freedom as equal). This implies 



A/h = AfH,cq(fcoqi?)' ( ^ 

5* 



1/3 



(24) 
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FIG. 1: PBH abundance as a function of horizon mass for 
power-law power spectra with n = 1.25 and n — 1.5 (left 
and right-hand sets of curves respectively) calculated using 
the Press-Schechter formalism with Ath = 1/3 (solid line) 
and the peaks formalism with (th ~ 0.7 and 1.2 (dotted and 
dashed respectively). 

In the early Universe g^, is expected to be of order 100, 
while 5+,oq ~ 3, fcoq = 0.07Slni/i^ Mpc~^. The horizon 
mass at matter-radiation equality is given by 

^H.oq - — 2prad,cq-Hcq = ^ , (^5) 

"'cq^cq 

where = 24000r2nift.^ and (assuming three species of 
massless neutrinos) rirad,o^^ = 4.17 x 10~^ so that 

MH,eq = 1.3 X lO^%Q^h^)-\ . (26) 

If we take n,^h^ = 0.14 then Mn^cq = 7 x lO^Og. 

In Fig. n we show various calculations of the abun- 
dance flpBii{M) for power-law primordial power spectra 
with spectral indices n = 1.25 and 1.5. The traditional 
calculation with Ath = 1/3 is compared with the peaks 
theory calculation for the two thresholds Cth = 0.7 and 
1.2. We see that the two peaks theory calculations actu- 
ally bracket the traditional calculation. The high value 
of Cth, corresponding to the lower abundance of PBHs, 
is the one which corresponds to peaks surrounded by a 
FRW Universe, and hence is likely to be more appropri- 
ate for the cosmological models under discussion. 

While we advocate use of the peaks theory expression 
Eq. H21|l to calculate the mass function, we see in the 
figure that the curves have similar shapes to those of the 
traditional calculation. In fact, if the peaks theory and 
Press-Schechter expressions were exactly the same, the 



thresholds would simply be related by Eq. (|13|l . which in 
radiation domination would give Ath = 4Cth/9. It turns 
out that this correspondence does hold quite accurately 
for our results even at the low abundances ripBH ~ 10~^° 
which are close to current observational bounds, breaking 
down only at much lower abundances where peaks the- 
ory is systematically higher than Press-Schechter. We 
therefore have quite a good correspondence: Cth = 1.2 is 
equivalent to Ath — 0.5, and Cth = 0.7 to Ath — 0.3. 



IV. DISCUSSION 

We have provided a new calculation of the abun- 
dance of PBHs generated by primordial density pertur- 
bations. By using a metric perturbation variable rather 
than the density contrast, a PBH formation criterion 
can be applied directly to the initial perturbation spec- 
trum. Within this formalism, we have found that the 
PBH mass spectrum is best computed using the theory 
of peaks, rather than the standard Press-Schechter-like 
calculation. 

Given the considerable uncertainties involved, our re- 
sults do not lead to any drastic revision of the PBH for- 
mation rate, but do put the calculation on a sounder 
theoretical footing. Our mass function can be fairly well 
approximated by that of the standard calculation in the 
region of interest (JIpbh ~ 10"^"), if the threshold den- 
sity Ath is taken in the range 0.3 to 0.5. This range of 
threshold values is however significantly lower than the 
value Ath — 0.7 suggested by the simulations of Niemeyer 
and Jedamzik |^, and in fact encompasses the value 
Ath = 1/3 used in the earliest PBH literature. How- 
ever, we advocate that anyone using our results adopts 
the peaks theory expression for the mass function given 
by Eq. 

After completion of this paper Ref. ^3 '^^^ brought 
to our attention. This paper uses the constraints on the 
metric perturbation variable ^ from Ref. J|] to calculate 
the PBH abundance, but does not use the peaks formal- 
ism. 
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